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The long-time and long-distance asymptotic behavior of the x spin correlations at finite tempera- 
ture in an anisotropic spin-1/2 XY chain is determined numerically. The decay of the correlations is 
exponential in both space and time. Similar exponential decay of correlations was already found ear- 
lier in the special case of the isotropic model, where analytical expressions for the decay rates could 
be derived via a mapping to a different model. While no such mapping is known for the anisotropic 
model, the asymptotic correlations can be very well approximated by a natural generalization of the 
known analytic results for the isotropic case. 
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Spin-1/2 chains are simple quantum many-particle sys- 
tems which can be defined in terms of a small number of 
coupling parameters but nevertheless offer a rich variety 
of interesting static and dynamic phenomena. The sim- 
ple structure of these models has made possible a large 
number of links to other models and fields and has fas- 
cinated researchers ever since the early days of Ising [l| 
and Bethe (H. 

The XY chain 0, Q] is an especially simple example 
since it can be mapped to a system of noninteracting 
lattice fermions. Its dynamics are nevertheless nontrivial 
since single-site spin operators are mapped to strings of 
Fermi operators and hence two-spin correlation functions 
correspond to many-fermion correlation functions which 
may become cumbersome to evaluate analytically. 

However, numerical calculations have been important 



in hinting at the direction where to look for exact analyti- 
cal results. An early example is the numerical calculation 
of Sur et al. Q for up to nine spins, suggesting a Gaus- 
sian behavior of the infinite temperature x autocorrela- 
tion of the isotropic XY chain. That numerical evidence 
was soon corroborated by independent rigorous proofs for 
the general XY model from two groups [6|-l8| . Here we re- 
port numerical results for the long-time and long-distance 
asymptotic behavior of the x spin pair correlation func- 
tions of anisotropic XY chains. Analytical results for 
these correlations were derived for the special case of the 
isotropic model only Q. Our numerical results indicate 
that the analytic formulae valid in the isotropic case pos- 
sess natural extensions into the anisotropic regime. 

The S* = 1/2 XY model 0, Q is defined by the Hamil- 
tonian 
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with anisotropy parameter < 7 < 1. The limiting 
cases 7 = and 7 = 1 are the isotropic XX and trans- 
verse Ising (TI) chains, respectively. The Jordan- Wigner 
transformation [sl. 13. [lo| 



S+ = <-''al , Sr = a,(-l)Er=i ^--^ ^ (3) 

between the spin-1/2 operators S-,Sf = Sf ± iSf at 
lattice sites i and the creation and annihilation operators 
a|,ai of lattice fermions maps the spin Hamiltonian ([1]) 



(2) to a Hamiltonian of noninteracting fermions: 



ajai+i + al^j^a, + 7 (alaj^^ + a, 
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Note that for 7 7^ the number of fermions is not con- served. From (l2]),([3|) it is evident that the spin corrc- 
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lation functions {S^{t)S^) are essentially fermion den- of the longitudinal correlation function (S'f(t)S'|^) is ^ t 



sity correlation functions, while the correlation functions 
{Sf{t)Sj) are much more complicated when expressed in 
the fermion representation. The fermionic identity 



(5) 



converts the string of signs appearing in ^ into an ex- 
pectation value of 2{i + j — 1) fermion operators, which 
may be expressed as a Pfaffian [ll[ by Wick's theorem. 
Pfafhans are close relatives of determinants and play a 
role in several statistical mechanical problems ■ Their 
numerical evaluation proceeds along similar lines as that 
of determinants. The numerical calculations of this study 
were all performed for spin chains with open boundary 
conditions. Periodic boundary conditions, while desir- 
able from an aesthetic point of view, lead to additional 
boundary terms in the fermionic model 0, 0] which make 
numerical calculations awkward, if not impossible. Of 
course, finite-size and boundary effects are a matter of 
concern in any numerical calculation. They are also a 
topic of research in their own right and have been stud- 



ied earlier Il3l Il4l |. In the present study, however. 



we focus on the asymptotic behavior of bulk spin corre- 
lations in the thermodynamic limit. To make sure that 
open-chain numerical results pertain to that situation, 
only spins sufficiently far from the boundaries of suffi- 
ciently long chains may be considered. We have checked 
that the numerical results to be presented below are not 
subject to finite-size or boundary effects during the time 
intervals shown. 

Several results about exact and asymptotic properties 
of the dynamic spin correlation functions have been ob- 
tained over the years. The asymptotic time dependence 
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for all i and j in the bulk of the system at all tempera- 
tures [i3-[i3|. This can be traced back to the properties 
of the one-particle density of states of the Jordan- Wigner 
fermions and the fact that {S^{t)Sj) is related to fermion 
density correlations. For i and j close to the boundary 
of a long open chain the situation is a little more compli- 
cated 
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Due to its more complicated structure in the 
fermion representation, the transverse correlation func- 
tion {Sf{t)Sj) is more sensitive to temperature varia- 
tions. At infinite temperature it vanishes for i ^ j and 
shows Gaussian decay for i = j At zero temper- 

ature that correlation function exhibits an asymptotic 
power-law decay in both space and time for the isotropic 
(XX) chain, jli, [H]. For finite temperature and in 
the isotropic case Its et al showed that the decay 
of {Sf{t)Sj) is asymptotically exponential in both space 
and time; numerical calculations could be used to as- 
sess the range of validity of the exponential asymptotics. 
Similar exponential behavior at finite temperature was 
also observed in the TI chain [23| and in rather general 
one-dimensional gapless integrable models [2lj. Asymp- 
totic finite-temperature correlations of the TI chain were 
also studied in (22l . [23| ; many interesting results on the 
zero-temperature correlations of that model were recently 
obtained by Perk and Au-Yang [23 |. 

Before discussing our numerical results it is useful to 
recapitulate what is known [25| about the ground state 
of the XY Hamiltonian ([T]) for general anisotropy 7. Em- 
ploying a Fourier transform followed by a Bogoljubov 
transform. H can be brought into diagonal frce-fermion 
form with the single-particle spectrum 



£/c = — sign(ft, + Jcos fc) Y (/i + Jcos fc)2 + 7^ sin fc; (— tt < fc < tt) 

I 



(6) 



The spectrum (for \h\ < \J\) has two branches with nega- 
tive and positive one-particle energies, respectively, with 
a gap of size Ae = 2|7|-\/ — at the critical wave 
vector given by cos fee = —j- The ground state in the 
fermionic picture has all negative-energy states occupied 
and all others empty. The spectral gap closes at the crit- 
ical field strength \h\ = \ J\ (for arbitrary 7) and for the 
isotropic case 7 = and arbitrary h (the XX chain), 
the case considered by Its et al Q . For h ^0 the ground 
state exhibits long-range order along the z axis; the mag- 
netization (5*1) is known analytically [2^. For 7 7^ and 
\h\ < \J\ the ground state also exhibits long-range order 



in the xy plane [2f 
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At finite temperature the equilibrium spin correlation 
functions {S^S"^^) were shown to decay exponentially 
with r for all a = x,y, z [2^. For the XX chain, 7 = 0, 
and subcritical fields, \h\ < |J|, the finite-temperature 
asymptotics of {Sf{t)Sf_^^{0)) were found to be exponen- 
tial in both space and time [9| . The derivation of that re- 
sult proceeded by mapping the correlation function to the 
solution of a classical nonlinear integrable system related 
to the nonlinear Schrodingcr equation. Despite that in- 
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tricate derivation the result can be simply interpreted Q 
in terms of the free energy of quasiparticlcs with disper- 
sion ^ (for 7 = 0). Generalizing to 7 7^ 0, we conjecture 
the following asymptotic formula: 



where 



exp (-/(?!, 0)) , n/vot > 1 
Cf^"'' exp(-/(n,i)) , n/vot < 1 



fin,t) = 



and where: 
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V = — In 
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tanh ■ 



while fco is determined by -^Ska ~ j- Furthermore, 

dsk 



vq = max 

fc|<7r 



dk 



is the maximum value of the group velocity at which the 
fermionic quasiparticles move. In the spacelike region, 
n > vot, the correlation is constant and equal to its sta- 
tionary value. This is due to the fact that the two spins at 
distance n have not been able to communicate with each 
other yet. Only later can the fermionic quasiparticles of 
energy Sk transfer any information at group velocity 
between the two spins. For zero magnetic field, h = 0, 
the maximum velocity vq turns out [27[ to be 



= (1 - 7)-^; 



(12) 



for nonzero field there is no simple closed expression for 
Vq, but numerical evaluation is simple. 

The crossover between spacelike and timelike regions 
is illustrated in Fig. [1] In that figure we show a loga- 
rithmic plot of \{Sfit)Sf^J\^ for n=4, 9, 14, 19 at tem- 
perature T/J = 2 and anisotropy 7 = 0.3. We observe 
that each function is almost perfectly constant up to the 
time Votn = n, where it bends smoothly into exponen- 
tial decay with weak superimposed oscillations. The de- 
cay time does not show any significant dependence on n. 
The inverse decay time predicted by ([H]) for the asymp- 
totic regime of the uppermost curve is given by the slope 
of the dot-dashed line and matches our data very well. 
The linear variation with n of the intercepts at t = 
in this logarithmic plot reflects the well-established ex- 
ponential decay of the equal-time correlation function 
(5f5f+„)~exp[-n/e(T)]. 

The inset to Fig. [T] shows again the curve n = 19 of 
the main plot along with curves for the same correlation 
function at different temperatures. Now the crossover 
between the space-like and the time-like regime occurs 
roughly at one common value of Jt. In the time-like 
regime, the slope changes from one curve to the next. 
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'^■'■^) FIG. 1: Main plot: Spin pair correlation function (Sf (t)5J) 
in the bulk regime of the anisotropic XY chain, 7 — 0.3, 
without external field, /i = 0, at T/J = 2, for j — 49 
and i = 30, 35, 40, and 45 (solid lines, bottom to top) in an 
open chain with A'^ = 100 spins. Plotted is the logarithm of 
(11) \{S^{t)Sj)\'^. The marks on the curves denote the boundary 
between space-like and time-like regions, compare (|8]). The 
inset shows the same function in the same representation, but 
here the sites are kept fixed {i = 30, j = 49) and the temper- 
ature is varied {T/J = 1, 10, 100, top to bottom). The slopes 
of the dot-dashed straight-line segments represent the decay 
rate inferred from the asymptotic expression (|8]). 



which reflects the T-dependence of the decay time, while 
the variable intercept in the space-like regime reflects the 
T-dependence of the correlation length. 

In Fig. [2] we show results for fixed spatial distance 
n — lA and fixed temperature T/J ~ 10, for varying 
anisotropy 7. For growing 7 the decay time in the asymp- 
totic regime grows, as does the crossover time between 
the space-like and time-like regimes. For 7 1 the 
crossover time diverges since the velocity vo ([T^ van- 
ishes. The decay time also diverges as 7 — >■ 1. This is to 
be expected, since in the absence of a magnetic field the 
model at 7 = 1 (i.e. the transverse Ising chain) does not 
display any dynamics. The 7-depcndence of the correla- 
tion length can be read off from the t = intercepts of 
the curves shown in Fig. [2] 

Fig. E] shows the autocorrelation function (n = 0) over 
a longer time interval, demonstrating very clearly how 
well the conjectured formula (|8]) fits the numerical re- 
sults. 

The temperature dependence of the autocorrelation 
function is displayed in Fig. |4] for fixed anisotropy, 
7 = 0.5. Two regimes can be clearly distinguished: at 
short times the autocorrelation is a Gaussian, crossing 
over to an exponential behavior at longer times. The 
crossover time seems to grow roughly logarithmically 
with temperature, so that at T — )■ 00 the well-known 
Gaussian behavior [B-Q emerges. 

The dependence of the spin correlation function on the 
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FIG. 2: Same as Fig. [T] for fixed distance n = 14 between the 
two spins, and fixed temperature, T / J = 10, aX h = Q. The 
anisotropy parameter is varied: 7 = 0,0.1,0.3,0.5,0.7 (bot- 
tom to top). The marks on the curves denote the boundary 
between space- like and time- like regions, compare ([S]). The 
slopes of the dot-dashed straight-line segments represent the 
decay rate inferred from the asymptotic expression 
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FIG. 3: Spin autocorrelation function {Si{t)Sf) in the bulk 
regime of the anisotropic XY chain, at /i = 0, T/J = 1, 
for i — 100 in an open chain with A'' = 200 spins. Plotted 
is the logarithm of \{Sf{t)Sf)f . The anisotropy values are 
(bottom to top) 7 = 0.1,0.3,0.5,0.7, and 0.9. The slopes of 
the dot-dashed straight-line segments represent the decay rate 
inferred from the asymptotic expression (|8]). 



external field was also observed to follow the asymptotic 
formula ([8]) for subcritical {\h\ < \J\) field values. 

To conclude, we have found convincing numerical ev- 
idence for exponential decay of the finite-temperature x 
spin correlations of anisotropic XY chains in both space 
and time. We suggest an analytic formula ([8]) for that 
exponential decay which fits the numerical data well and 
which generalizes the known analytical result 0] based on 
a mapping to a classical nonlinear intcgrable system. At 
high temperature our numerical results connect well to 
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FIG. 4: Same as Fig. [3] for fixed anisotropy 7 = 0.5 and 
varying temperature T = 1,10,100,1000,10*, and lO'^ (top 
to bottom). 



the known [6|-|8| Gaussian infinite-temperature behavior. 
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